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Abstract: We compute the NLO QCD corrections to the pair production of VF-bosons in 
association with two jets at the Tevatron and the LHC. This process is an important back- 
ground to heavy Higgs-boson production in association with two jets, either in gluon or 
weak boson fusion. We consider leptonic decays of W-bosons and include all the spin corre- 
lations exactly. For natural choices of the renormalization scale, the NLO QCD corrections 
to pp(p) — >• W + W~jj are moderate but different for various values of the center-of-mass 
collision energy at the LHC and the Tevatron, emphasizing the need to compute them 
explicitly. 



Keywords: 



NLO Computations, QCD, Jets, Hadronic Colliders 



Contents 



[lj. Introduction |] 

§. Technical Details | 

|2.1| Processes with W + W~ pair, a quark pair and two gluons || 

2.2 Processes with W + W~ and two quark pairs [?] 



|2.3| Checks on the calculation 

^ Phenomenology 11 
|3.1| Results for the Tevatron 11 



3J Results for the LHC |14 
^j. Conclusions 17 
|A]. Results at a fixed phase space point 20 



1. Introduction 

The Large Hadron Collider (LHC) has begun to explore the Standard Model (SM) of 
particle physics in a new energy regime and will in time gather more data than any previous 
hadron collider experiment. To further our understanding of the SM and as to what may 
lie beyond it, we attempt to describe outcomes of proton collisions in sufficient detail, for 
comparison with the observed data. An accurate knowledge of SM processes is particularly 
important as their cross-sections are often much larger than those for many interesting 
New Physics processes. Unless physics beyond the SM presents itself in a stark way, 
disentangling it from SM backgrounds will require an accurate description of the latter. 
Parton level calculations at leading order (LO) in the strong coupling constant are often 
insufficient for this purpose. They exhibit a strong unphysical dependence on factorization 
and renormalization scales, leading to large uncertainties in the predictions. Data-driven 
estimates of the backgrounds are also subject to large uncertainties if they rely on LO 
theoretical predictions: here the idea is to determine the normalization of the LO cross- 
section for a given background process in a region essentially free from any New Physics 
signal. Once the LO is "validated" using data, one extrapolates it to the region of interest. 
It is clear that such a procedure can only work if higher order QCD corrections are uniform 
over phase space, which is not guaranteed in general. As follows from many successful 
analyses at the Tevatron, a good way to reduce the uncertainty is to extend the theoretical 
description of a given process to next-to-leading order (NLO) in perturbative QCD. 

The past five years have seen an extraordinary progress in the development of meth- 
ods that are suitable to deal with NLO QCD computations for high-multiplicity processes. 
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Refinements of traditional computational techniques based on the Passarino-Veltman re- 
duction of tensor integrals led to the development of highly-efficient, Feynman-diagram- 
based technology for NLO QCD computations |lj, ^ ||||. At the same time, new techniques 
based on unitarity and on-shell methods |5|, q, |7], ||, 10 1 sufficiently matured to be- 
come relevant for practical applications. As a result, a large number of 2 — > 4 processes 
were studied at NLO in QCD in the past two years. The list includes pp — > W(Z,j) + 3 
jets fll], 0, |§ 0, |15|, [ft 



pp 



ttbb 17, 



pp 



tt + 2 jets 



bbbb m 



pp^tt -)• (2|, H, and -)■ W + W + + 2 jets [|J. This last process has been 

implemented recently in POWHEG-BOX [Pq| . This combines NLO accuracy with a parton 
shower detailed description of the final state. 

The first 2 — > 5 process, pp — > W+4 jets, has also been computed recently through NLO 
QCD using on-shell methods (2(|. Some groups also started employing those advances with 
the intention of developing a platform for fully automated NLO QCD calculation [27, 28]. 

The goal of this paper is to study the production of a pair of TV-bosons in association 
with two jets in hadron collisions, including the NLO QCD corrections. The production 
of a VF-boson pair in association with zero, one or two jets is an important background 
to searches for intermediate and heavy Higgs boson, where the decay H — > W + W~ opens 
up. At the Tevatron, searches for intermediate-mass Higgs bosons treat processes pp — > 
H + n jets, n = 0, 1, > 2 separately, because dominant backgrounds depend on the number 
of identified jets in the final state (see e.g. |2£|). While most of the sensitivity in Higgs- 
boson searches comes from the process with the largest cross-section, pp — > H + jets, the 
production of the Higgs boson in association with two jets is also relevant |3^, plj. Because 
pp — > W + W~ jj is an irreducible background to the Higgs-boson production in association 
with two jets, it is important to have NLO QCD predictions for this process. 

There is yet another reason to want an improved description of W + W~jj production 
in hadron collisions. At the LHC the Higgs boson can be produced with a sizable cross- 
section in weak boson fusion (WBF) |30|, In addition to the Higgs-boson decay products 
which, as we assume, are pairs of VF-bosons, the signature of the process involves two 
forward tagging jets. In this case, pp —> W + W~jj is the irreducible background. The 
Higgs-production cross-section in WBF is known through NLO QCD [| 



p3| , |34j, and 

it is desirable for the dominant background process to be known to the same order in 
perturbative QCD as well. 

Finally, we note that jets, charged leptons and missing energy is one of the classic sig- 
natures of dark- matter- type processes at colliders. In such scenarios the missing energy ap- 
pears due to the dark-matter candidate escaping the detector. The process pp — > W + W~jj 
is a SM background with a similar signature, where leptonic decays of VF-bosons lead to 
invisible neutrinos. 

Several studies in the past decade addressed the production of VF-boson pairs in hadron 
collisions, including NLO QCD corrections. In particular, pp — > W + W~ with no jets was 
studied in Refs. [35, 36, 37]. The production of a pair of H^-bosons in association with 
one jet including decays to leptons was studied through NLO QCD in Refs. [38, p9| . In 
both cases, for the choice of the renormalization and factorization scales \x = M\y, QCD 
corrections were found to be significant, of the order of (25 — 50)%. These results further 
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motivate the need to understand the production of W + W~ in association with two or even 
more jets at NLO in QCD. 

In this paper we allow for leptonic decays of the TV-pair including all spin correla- 
tions. Dilepton final states are the ones that are relevant for ongoing Higgs searches at 
the Tevatron and, in general, these final states provide the cleanest signature to identify 
the production of TV-bosons at a hadron collider. For this reason, we find it reasonable to 
focus on these states only. 

We remind the reader that the branching fraction for the TV-boson to decay to a 
definite-flavor lepton final state is about ten percent. Since we have two TV-bosons decaying 
leptonically, we get a hit by a factor 0(1O~ 2 ) when the cross-section for dilepton final state 
is compared with the cross-section for stable TV-bosons. It is therefore amazing that the 
cross-section for the process pp — > (TV + — > f^fJ> + ) + (TV~ — > e~v e ) + jj is still reasonably 
large. In particular, we find that the cross section for the LHC running at an energy of 7 
TeV is around 40 fb, which means that a few of these events should have already been seen 
at this collider at the time of publication, and quite a significant number of such events 
should be produced at the LHC by the end of the next year. The cross section further 
increases to about 0.14 pb at 14 TeV, so there is no doubt that the experimental study 
of this process is feasible. Even at the Tevatron, where the cross-section with the "Higgs- 
like" cuts for one flavor assignment is just 2.0 fb, assuming fifty percent efficiency, about 
40 e + e~vvjj, fi + pTvvjj , fi + e~~vvjj, e + [i~vvjj events should have been recorded already. 

The computation of NLO QCD corrections to hadro-production of W + W~ , in associ- 
ation with two jets, is also interesting from the point of view of further developing on-shell 
methods for one-loop computations. Recall that, as currently formulated, on-shell meth- 



ods require ordering of external lines which is achieved by working with color-ordered [40] 



or primitive amplitudes [41, |42j, 43 1. These techniques work best if all external particles 
carry color charges, while their implementation becomes more involved as the number of 
colorless particles in the process increases. The only process which involves two colorless 
particles that has been computed with unitarity methods before, pp — > W + W + jj p4|j . 
is simpler than the calculation presented here since, among other things, only a smaller 
number of sub-processes contribute. As explained in the next Section, the presence of two 
colorless particles in the final state whose total electric charge is zero provides some ad- 
ditional difficulty. Nevertheless, it is possible to handle these complications with on-shell 
methods. 

The remainder of the paper is organized as follows. In Section B we provide technical 
details of the calculation. In Section |3] we discuss phenomenological results for the QCD 
production of W + W~jj, with leptonic decays of the TV-bosons, at the Tevatron and the 
LHC. We conclude in Section 4. We provide numerical results for various one- loop primitive 
amplitudes, as well as squared amplitudes summed over helicity and color for (TV + — > 
z/ At /x + ) + (TV - — > e~u e )jj hadronic production in the Appendix. 

2. Technical Details 

In this Section we present technical details specific to this calculation. Within a subtraction 
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formalism, a NLO calculation involves three components: virtual corrections, real emission 
corrections, and subtraction terms for soft and collinear divergences. The virtual correc- 
tions are computed using D-dimensional generalized unitarity & Pj. A detailed description 
of the implementation of this method can be found in [44]; we have followed this imple- 
mentation, modifying and extending it to deal with the presence of an additional P^-boson 
in the final state. 

The full one-loop amplitude can be built by summing products of color-ordered partial 
amplitudes (il], 42, 43 1 over all permutations of the colored particles, with appropriate color 
factors. The partial amplitudes are further decomposed into primitive amplitudes. The 
ordering of all particles with color charges is fixed in primitive amplitudes. D-dimensional 
unitarity cuts reduce one-loop primitive amplitudes to linear combinations of products of 
tree-level helicity amplitudes, which are computed using Berends-Giele recursion relations 
1 45]. These relations are also used to compute tree- level amplitudes which are required for 
calculations of LO cross-sections, real emission corrections and subtraction terms for soft 
and collinear emissions. We implement subtraction terms following the Catani-Seymour 



procedure [j4(|, with the a-parameter optimization as described in Ref. [47, 48]. We embed 
our calculations within the framework of the MCFM program [49] and use the QCDIoop 
program to calculate the scalar one-loop integrals 0] . 

Since only color-charged particles are ordered in primitive amplitudes, all possible 
insertions of the Vy-bosons must be considered when tree-level or one-loop primitive am- 
plitudes are computed. While this implies a certain amount of non-trivial book-keeping in 
the construction of a numerical program, this can be done without much trouble. The real 
problem, however, is that cuts of different parent diagrams must be combined in certain 
cases to produce gauge-invariant tree-level amplitudes in the context of unitarity cuts. This 
implies that different parent diagrams can not be treated independently and this creates 
considerable overhead. Furthermore, we must include the possibility of the W + W~ pair 
being produced via an intermediate neutral vector boson, such as an off-shell 7 or Z. 

In this calculation, we do not consider the production of top quarks in the final state, 
as these are processes with a distinct experimental signature. Furthermore, we neglect top- 
quark contributions in virtual diagrams and treat all other quarks as massless. Since top 
quarks in virtual diagrams originate from b — > Wt* transitions, we decided to completely 
exclude bottom quarks in our calculation as well. This is a reasonable approximation 
since the 6-content of the proton is subdominant both at the Tevatron and the LHC. 
We also neglect g* — > bb splitting, both real and virtual. We believe that this effect is 
also quite small as can be seen from the 6-quark contribution to the QCD /3-function, 
relative to contributions of gluons and four other quarks. Although we do not expect 
that the complete omission of quarks in the third generation impacts our results in any 
significant way, we hope to include them in the calculation in the future. The framework 
to do so, within the generalized D-dimensional unitarity approach, has already been fully 
elaborated in Refs. [51, 52]. Before continuing, we point out that in this paper we do not 
include contributions from one- loop diagrams where the j/Z or the W + W~ pair couple 
directly to a loop of virtual quarks, creating a diagram of the "light-by-light scattering" 
type. These diagrams form a finite, gauge-invariant class of amplitudes that can be dealt 
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with separately. In particular, amplitudes for the partonic process gg — > W + W gg which 
does not appear at tree-level also belong to that class of amplitudes. As pointed out in 



Ref. [53], processes of that type may be quite important because of the large gluon flux at 
the LHC. We plan to return to the discussion of the amplitudes where pair couples 

directly to a closed quark loop in a separate publication. Finally, in this calculation we 
neglect mixing between up and down quarks of different generations and set the CKM 
matrix to the identity matrix. 

The production of W + W~jj can occur through both electroweak and QCD mech- 
anisms. The NLO QCD corrections for the electroweak production have already been 



calculated in Ref. [54|. While these mechanisms can interfere even at leading order, these 
interference terms are strongly suppressed. First, at partonic level, the electroweak pro- 
duction of W + W~jj involves four quarks. However, given the large gluon luminosity at 
the LHC, four-quark contributions to W + W~jj production cross-section amount to only 
about 15%. Moreover, the interference that occurs in a four-quark process can only happen 
for certain combinations of quark flavors, and it is color-suppressed. We therefore neglect 
this interference, and present results for the QCD production alone. 

In order to describe tree-level and one-loop virtual corrections to pp(pp) — > W + W~jj 
we require partonic processes with either two quarks and two gluons — > q\q 2 gg W + W~ or 
with four quarks — > q±q 2 q^qA W + W~. Given the difference in color- and flavor structures, 
we discuss these two partonic processes separately in the next two Subsections. 

2.1 Processes with W + W~ pair, a quark pair and two gluons 

In this Section, we consider the partonic process — > q~iq 2 ggW + W~~ . Since we neglect 
mixing between up and down quarks of different generations, in two-quark amplitudes 
both quarks q± i2 have the same flavor. 

To obtain the full NLO cross-section, we need to consider all possible crossings between 
the partons; the initial state partons - as well as the jets - may be either gluons or quarks. 
The tree level amplitude for the process — > (qq) + (W + — > v fl +n + ) + (W~ — > e~ +v e )+g+g 
can be written as 



-4 tree (<7i,<72; ,zv/i + ,e ,v e ;g3,94) = gl p w(s u ^+)Pw(s e -vJ 

X ((T a:i T a ^- lli2 A (q 1 ,q 2 ;g 3 ,g 4 ) + {T a *T a *) lxi2 A Q {q x , q 2 ; g 4 , g 3 



(2.1) 



In Eq.(fO|), g s and gw are the strong and weak coupling constants, respectively, leptonic 
labels have been suppressed on the right hand side, and Pw are Breit-Wigner propagators 

Pw(s) = M2 * rM - (2-2) 

s — Myy + iTwMw 

with siu = (pi + Pv) 2 ■ In addition, My/ and T\y are the H^-boson mass and width, and 
the generators of the SU(3) color group are normalised to Tr(T a T b ) = 5 ab . In Eq.(fO]), 
Aq denote the color-ordered amplitudes. The flavor of the quark line fixes the electric 
charges of q\ and q 2 and, simultaneously, the ordering of W + and W~ along the quark 
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Figure 1: Primitive amplitudes a) A 1 (q 1 ,q 2 ,g 3 ,g4), b) A 1 (q 1 ,g 3 ,q 2 ,g 4 ), c) Ai(gi, # 3 , g 4 , q 2 ), and 
d) A^ 2 ' (51,(72,53,34)- W bosons are not shown. 



line. However, as we pointed out before, the relative ordering of W bosons and gluons is 
not fixed. Additionally, we need to consider the possibility that VF-bosons are produced 
through an intermediate (off-shell) Z-boson or photon. Thus we write 

4)(gi,<&;^,s4)=4 w) (&>^^ (2.3) 

where the first term describes an amplitude where H^-bosons couple directly to the quark 
line and the second term describes an amplitude where such coupling occurs through a 7 
or Z. The factor c( q2 ' h ^ is given by 

C ( q ,h) = 2Q (q) gin 2 Qw + p z ( Sz )(tW _ 2Q(«) S in 2 Ow) , (2.4) 

where and h are the electromagnetic charge and helicity of the quark q, T3 ^ = 1 and 
T 3 (+) = 0, 6w is the weak mixing angle, and sz = {pw+ +Pw~) 2 = (Pi> +Pfi+ +P e - +Pv e ) 2 - 
Note that, because VF-bosons only couple to left-handed quarks, the first term in Eq. fl2.3| ) 
is zero if the quark is right-handed. We account for the decay — > l^{pi) + viijpv), by 
using the polarization vectors constructed from lepton spinors. For example, in case 
of the W + boson, the polarization vector reads 

AW \ U{p v )^^-V{pi+) I-75 

e-{Pu,Pi+) = — -, — 72 — ; 7- = n • (2-5) 

[Pi+ + PvY 2 

The computation of real emission contributions requires tree amplitudes with an ad- 
ditional gluon in the final state. The color decomposition reads 

A tICC (qi,q2;v ll ,n + ,e- ,v e ;g 3 ,g A ,g 5 ) = g^ ( ) Pw(s v „+)Pw(,s e - Ve ) 

W2J (2.6) 
x £ (T a ^T a ^T a "sy iii2 A (q 1 ,q 2 ;g (T3 ,g a4 ,g (T& ), 

where Si denotes the permutation of i indices. The flavor /helicity properties of the ampli- 
tudes with three and two gluons are identical and have been discussed for the two-gluon 
case. 
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The decomposition of the one-loop amplitudes in terms of left-handed primitive am- 
plitudes [gl], Q reads 

A 1L {qi,q 2 ;^,fJ- + ,e~,iy e ;g3,g4) = 5s( -7= J p w(s u „+)P w {s e -pJ 



V2, 

+ (T^T-3T^) ?i . 2 (r-) :ci ^^ 1 (g- 1 ,^3^2,^) (2 ' 7) 
+ (^ 2 ^ 1 )H J2 (r CT3 r CT4 ), 1 , 2 ^(«i^ 2 ,^ 4 ,9 1 



X 

<TGS 2 



0-3 / 

+ ^-Gr 4 ^^ 1/2l (gi,g 2 ,5' CT 3'5a 

In Eq.( |2.7| ) we introduced the color factor 

Gr 4 = N c (T a ^T a ^)- i . - Ti{T a ^T a "i)5- ili2 . (2.8) 

We build up the virtual amplitude from eight primitive amplitudes: A\{q\, 54,53,92), 
(91,53,92,54), ^1(91,92,54,53) and ^' 1/2] (9i,92,53,54), shown in Fig. [j] and another four 
amplitudes, obtained by swapping the gluons 53 o 54. In Fig. [j], we introduce a 'dummy 
line' for the primitive amplitude A^ 2 ' (qi, q 2 , 53, 54). This allows us to draw this primitive 
amplitude - which has the external gluons attached to a fermion loop - as formally having 
six loop-momentum-dependent propag ators 1 . The W-b osons couple to the dummy lines, 
but dummy lines cannot be cut. 

2.2 Processes with W + W~ and two quark pairs 

We now consider the case of amplitudes involving two qq pairs and the W + W~ pair. We 
first discuss the color and flavor structure of the tree-level amplitude — > (91929394) + 
(W + — >■ Vu + n + ) + (W~ —7- e~ + v e ) treating all particles as being in the final state. This 
process is described by Feynman diagrams with two continuous fermion lines connected 
by a gluon exchange, with VF-bosons being emitted from either of the two quark lines. 
Depending on the quark flavors and on the way the VF-boson emissions occur, we may 
have to assign quark fields in two different ways to the fermion lines: [9192], [9394] and 
[9194], [9392]- We refer to the first assignment as the "s-channel amplitude" and to the 
second assignment as the "i-channel amplitude" , see Fig. [| 

We begin by considering the s-channel tree-level amplitude. In this case, the color 
decomposition reads 

£ tree (9i,92;93,94;zv/^ + ,e~,^ e ) =5s ( ^7= ) Pw{s e+Ve )P w (s^- u ) 

W2J ' (2.9) 

x {kxuksii ~ Jfk^isiJ^iQU 92! 93,94)- 

We can further split the Bq amplitude into two separate types. The amplitude of the first 
type appears if a quark line radiates both T^-bosons and the other quark line radiates none. 



1 We count external VK-bosons, which are not shown in Fig. |l| 
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93 




(s-type) (t-type) 

Figure 2: Amplitudes of s-type (left) and of t-type (right), for the partonic process — > 
W + W~ gi 329394- The VF-bosons are not shown. 



The VF-boson can be radiated either directly from the quark line, or through an exchange 
of an intermediate ~f/Z. The amplitude of the second type arises when one W-boson is 
radiated from each of the quark lines. Examples of the corresponding contributions are 
shown in Figs. |3],[| for specific flavor assignments; it is clear how this classification generalizes 
to other flavors. 

We begin by discussing amplitudes of the first type. Since we set the CKM matrix equal 
to the identity matrix, flavors of fermions can not change along the fermion lines, so that 
flavors of q\ and q 2 as well as of q 3 and (74 are equal. Thus, for a set of flavor assignments, for 
which this contribution is allowed, there are four diagrams that contribute to the amplitude 
Bq. Examples are shown in Fig. ||[ We write the color-ordered amplitude as 

A)(9i,92;93,94) = B^ VW) {[q l ,W,W,q 2 },[q 3 ,q 4 }) + B { ^ W) {[q u q 2 },[q 3 ,W,W,q A ]) 
+ C^B^ /Z \[q ul /Z, q 2 ], [q 3 , q 4 }) + B^ /Z) ([q 1} q 2 ], [q 3 , 7 /Z, q 4 }) , 

where h 2 ^ = {— ,+} are the helicities of quarks q 2 and q^. We note that since VF-bosons 
couple only to left-handed quarks, the first term in Eq.( 2.1C| ) is zero for h 2 = 1, and the 



second term is zero for J14 = 1. The factors C^ q,fl ' are given in Eq. ( |2.4[) . The second term 
in Eq. ( |2.10| ) can be obtained from the first term by swapping momenta pq 1 o pq z and 
p q2 «->■ p q4 . The same swap can be used to obtain the fourth term in Eq. ( |2.10| ) from the 
third. 

We turn to the discussion of the amplitudes of the second type, which correspond to 
the emission of the W + boson off one quark line and the W~ boson off the other quark 
line. As a result of the emission, flavors change along each fermion line. An example of a 
diagram contributing to this amplitude is shown in Fig. |4|. As there is no contribution of 
the neutral vector boson in this case, the amplitude is nonzero only for h 2 = h.4 = — 1. The 
choice of flavors for ^1,^3 determine which VF-boson is radiated from which quark line. 

According to the flavors of the four quarks, only one of the s- or t-channel amplitudes 
can contribute, or both. Since the ^-channel amplitude is obtained by replacing q 2 q^ 
in the s-channel amplitude, everything that has been said about the latter applies to the 
former. Note that the replacement q 2 (74 also involves color indices, so that non-trivial 
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Figure 3: Sample tree-level diagrams for Bq(u,u, c, c). When both W^-bosons couple directly to 
the quark line, the flavors of the quarks determine the ordering of the VK-bosons. 



color-correlations appear in the interference of s- and i-channel amplitudes when both are 
allowed by flavor. 

For the computation of real emission corrections we need four-quark amplitudes with 
additional gluon in the final state — > (91929394) + (W + — y + fi + ) + (W~ — > e~ + v e ) + g. 
It is clear that the presence of an additional gluon does not modify the separation of 
amplitudes into s- and f-channel amplitudes, so that much of what has been said about 
the tree-level amplitudes remains applicable. In particular, the flavor structure is identical 
to the tree-level case discussed above. On the other hand, the color decomposition differs. 
For instance, for the s-channel amplitude, it reads 

4 



i3 tree (<?i,g 2 ,<?3,<?4,c/;^,M + ,e ,v e ) = gl(^^j Pw(s u ^+)P w (s e - De 
S i 3 i2 T f 1 i 4 B o(Qi, ( l2,q3,q4,9) + j^5i 3i4: T( ii2 Bo(qi,g,q2,q3,q4 



+ Si 1 i i Tf +3i2 B (q 1 ,q 2 ,g,q3,q4:) + jf 5 hi 2 T t 3 u B o{qu 92, 93, 9, 94) 



(2.11) 



Similar considerations apply to virtual corrections but the color decomposition is more 
involved in this case. For the s-channel virtual QCD amplitude it reads 

B 1L (qi, 92, 93,94! ZV^ + , e ~> p e) =5a(-7=) P w(s u ^+)P W (s e -pJ 

\V2/ (2.12) 
x {kwA^Bl (q^ ^2, 93, 94) + $i li2 Si 3 i 4 B[ 2) (91,92,93, 94)) • 
The amplitudes in Eq,( |2.12j ) are written through primitive amplitudes as 



Bi ] (91, 92, 93, 94) = (N c - -^-)B\ a '(q 1 ,q 2 ,q3,q4) - -^-B ( 1 a) (q 1 , q 2 , 93, 94 



>(«) 



(«), 



N, 



-^<(9- 1 ,92,93,94) 



B[ c) (91 , 92, 93, 94) + n f B\ a > (91 , 92, 93, 94 



(<*), 



(2.13) 



N, 
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Figure 4: Sample tree-level diagram for B (u,d,s,c). 



and 

B { i ] (<?i,<?2,<73, Qa) = (91,92,93, 94) + (l + ^)£i a) (9l, 92, 93, 94) 

1 C 1 n (2 - 14) 

+ jpBf ) (q l ,q 2 ,qz,qA) + j^, B x (9l, 92, 93, Qa) ~ jf B i (?l> 92, 93, 

Parent diagrams for the primitive amplitudes B^ ,b,c,d ' are shown in Fig. [|. The only 
primitive amplitude that receives contributions from six-point one- loop diagrams is B^ a \ 
Primitives and are simply the Born amplitudes dressed by a gluon loop on one of 
the quark lines, while fl" corresponds to a fermion loop contribution. For convenience, 
we again use dummy lines in Fig. ^; they allow us to consider every primitive amplitude as 
having a parent diagram with (formally) six propagators. We recall that VF-bosons couple 
to dummy lines, but that these lines cannot be cut. Berends-Giele recursion relations are 
modified in these cases to ensure that the correct primitive amplitudes are recovered. 

2.3 Checks on the calculation 

Various checks were carried out at all stages of the calculation. The squared matrix el- 
ements for the leading order and real emission processes were checked against Mad Graph 



1 55 ] for a few phase space points. This was done for all flavor combinations and all initial 
state parton configurations. Gauge invariance of various amplitudes was checked for both 
the external gluons and the VF-bosons (artificially setting the masses of the latter to zero), 
at leading and next-to-leading order. The subtraction terms of the Catani-Seymour dipole 
method were checked to cancel with the real emission terms in the limit when emitted 
partons become soft and/or collinear. We checked the double and single infrared poles of 
the virtual contribution, both at the level of primitive amplitudes and at the level of vir- 
tual matrix elements squared. These terms were also checked to cancel with the integrated 
dipoles. We also checked the independence of the cross-section of the a-parameter [47, 48]. 



Finally, the full one-loop amplitude is checked against an OPP-based, but otherwise com- 
pletely independent diagrammatic computation, at a few phase space points. We note that 
over six hundred Feynman diagrams are involved in a such a calculation. 

By default, our calculation was performed in double precision. For each phase space 
point, the double and single poles were checked against the analytically known results, and 
the coefficients of the OPP expansion were checked to have accurately solved the system 
of linear equations. If either of those checks failed, the amplitude at that phase space 
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Figure 5: Parent diagrams for one-loop primitive amplitudes B\ 



(a,6,c,rf) 



for (gig^^+W+W - 



where the flavors of the quarks are not specified. The VF-bosons are not shown. Shaded areas 
represent dummy lines which are not cut. 



point was recalculated using quadruple precision. We found that around 0.4% of primitive 
amplitudes had to be recalculated this way. 



3. Phenomenology 

In this Section, we discuss phenomenological aspects of W + W~jj production at the Teva- 
tron and the LHC. At the Tevatron, this process is a background to Higgs-boson production 
in association with two jets. We employ set of cuts discussed in the context of the Higgs- 



boson search in Ref. [29| and study related phenomenology. At the LHC, we consider 
the collision energy of 7 TeV and we show that the number of dilepton events related to 
W + W~jj production is sufficiently large to study this process in detail. 

Before moving on to a dedicated discussion, we briefly describe general features of our 
computation. The W-bosons are always produced on mass-shell and decay leptonically 
VF + M^ _ — > VnH z~v e . We note that, neglecting non-resonant contributions, the results for 
all lepton flavors = {e + e~, e + //~, /i + e~, n + fi~} can be obtained by multiplying our 
results by four. 

The mass and width of the VF-boson are taken to be M\y = 80.419 GeV and T\y = 
2.141 GeV, respectively. The width of the Z-boson is taken to be Y z = 2.49 GeV. The prop- 
agators for these particles take the Breit-Wigner form. The electroweak gauge couplings 
are computed using a QE D(M z ) = 1/128.802 and sin 2 6 W = 0.2222. We use MSTW08LO 
parton distribution functions for leading order and MSTW08NLO for next-to-leading or- 
der computations |p6| . The strong coupling constant a s (Mz) is part of the MSTW fit. It 
equals to 0.13939 (0.12018) at leading- and next-to-leading order, respectively. 

3.1 Results for the Tevatron 

By the end of Run II, the Tevatron will have collected just over 10 fb _1 of data for use 
in the search for the Higgs boson. At the very least, the two Tevatron experiments will 
be able to improve upon the exclusion limits for the Higgs bosons presented earlier in 



Ref. [57]. The search strategy is to separate relevant processes, depending on the number 



of jets produced with the Higgs boson. As follows from the analysis in Ref. [58], ten percent 
of all events with the Higgs boson at the Tevatron contain two or more jets. The process 
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Figure 6: The dependence on renormalization and factorization scale of the cross section for 
pp — > u^u + e ~v e jj at i/s = 1.96 TeV, where p — fiR = [if- Predictions at both LO and NLO in 
QCD are shown. 



pp — > W + W~jj is a SM background of significant importance. In Ref. |3g], the NLO QCD 
cross section for Higgs +2 jet production with the decay H — > W~(— > fi~r , fJi )W + (-^- z^ e e + ) 
is calculated. For the Higgs-boson mass of 160 GeV, the cross-section value ctnlo = 0.2 fb 
is found (we do not show the uncertainties which are significant). This cross-section is 
obtained with cuts that are similar to those used by the CDF collaboration in their Higgs- 
boson search |29[] . Specifically, jets are defined using the fej_-algorithm, with ARj 1 j 2 > 0.4. 
Jets must have p±j > 15 GeV and must be in the central region of the detector, \rjj\ < 2.5. 
It is required that two leptons, one with transverse momentum p_\_^ x >20 GeV and rapidity 
|% | < 0.8 and the other with transverse momentum p± t i 2 > 10 GeV and rapidity \rji 2 \ < 1.1, 
appear in the event. The invariant mass of the lepton pair is required to be larger than 
m hh > 16 GeV. Both leptons must be isolated. The specific requirement to this effect 
is that any jet within AR = 0.4 of a lepton must have a transverse momentum which is 
smaller than 0.1 p±j. The CDF collaboration uses a particular constraint on the missing 
transverse momentum. They introduce a function 

rZ CC =^sin[min(A<A,|)], 

with A<j) being the angle between the missing transverse momentum vector and the 
nearest lepton or jet. An event is accepted if ec > 25 GeV. 

We present NLO QCD results for the process pp — > W + (— > u fl fi + )W~ (— > e~v e )jj, at 
-y/s = 1.96 TeV, using the kinematic cuts that we just described. This allows us to study 
this process as a background to the Higgs-boson production. In Fig. ^ we show the scale 
dependence of the cross section for the process pp — > W + (— > W~{— > e~i> e )jj, both 
at LO and NLO in perturbative QCD (pQCD), with the scale ranging between M\y/2 and 
2M W . 

The leading order cross-section is ctlo = 2.5 ± 0.9 fb. This result is interesting since 
its uncertainty alone exceeds the cross-section for the production of the Higgs boson in as- 
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Figure 7: Kinematic distributions showing the opening angle between the leptons, 4> e - ti +, and the 
difference in rapidity of the two hardest jets, for the process pp — > W + (-^ v^ii + )W~ (-^? e~v e )jj at 
the Tevatron running at t/s — 1.96 TeV. The bands show renormalization and factorization scale 
uncertainty for Mw/2 < M < 2Myy, and the solid line is the prediction for u = M\y- 



sociation with two jets by about a factor between four and five. Clearly, there is no way to 
discuss observation of the Higgs boson in this channel unless the theoretical uncertainty 
on W + W~jj is improved. The situation indeed improves once NLO QCD corrections are 
computed. We find ctnlo = 2.0 ±0.1 fb - a significant reduction in scale uncertainty. How- 
ever, even after that reduction, we find that the uncertainty on the W + W~jj production 
cross-section is very much comparable to the absolute value of the Higgs-boson production 
cross-section in association with two jets. For this set of cuts, the NLO QCD computations 
lead to a prediction of about 80 efi + jj, e + fijj, e + e~jj, fi + fT jj events during Run II, using 
the discussed set of cuts and assuming 100% efficiency. 

There are other kinematic variables that one can use to improve upon a discrimination 
between the Higgs-boson production and the W + W~ production. For example, the opening 
angle of the two leptons is of particular interest. Indeed, if a pair of VF-bosons is produced in 
the decay of a scalar particle, their spins are anti-correlated. As a result, leptons from their 
decay tend to have small relative angles in the transverse plane. The ^> e -^ distribution in 
the case of QCD W + W~jj production is shown in Fig. |7|and the leptons are seen to have a 
preference to be back-to-back, in strong contrast to the Higgs-boson signal. No noticeable 
shape changes occur when the QCD corrections are included. In the second pane of Fig. [7|, 
we plot the rapidity difference between the two hardest jets Ar/j U - 2 = rjj 1 — r]j 2 , which is 
peaked at zero and falls off rapidly, with an almost vanishing fraction of the cross section 
having a magnitude of rapidity difference greater than four. Note that a requirement 
IAj^JjI > 4 is imposed when the Higgs boson is searched for in weak boson fusion. 

Finally, Fig. || shows the transverse momentum distribution of the charged lepton and 
Ht,tot defined as the scalar sum of the transverse momenta of all visible particles present 
in the final state plus the missing transverse momentum, -Ht,tot = J2j P±,j + P±,fi+ + 
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Figure 8: Kinematic distributions showing the transverse momentum of a lepton and Ht,tot, for 
the process pp —> W + (—> v^n + )W~ '(— > e~v e )jj at the Tevatron running at y/s = 1.96 TeV. The 
bands show renormalization and factorization scale uncertainty for Myy/2 < [i < 2Mw, and the 
solid line is the prediction for fi — M\y- 



p± e - + p_i_,miss- It follows from Fig. |8| that the shape of lepton transverse momentum 
distribution does not change but the -Ht,tot distribution becomes somewhat softer at 
NLO QCD. 



3.2 Results for the LHC 

The LHC is set to run at 7 TeV until the end of 2012, collecting 2-5 fb" 1 of data. As 
a result, a non-negligible number of dilepton events, originating from W + W~jj, will be 
observed at the LHC during this and next year, which warrants a phenomenological study 
of this process. The importance of pp — > W + W~jj process as a background to Higgs- 
boson production has been discussed extensively in the literature (see e.g. [^9|), especially 
with reference to the weak boson fusion production mechanism, where designed cuts on 
the jets can dramatically boost the signal to background ratio. In this paper, we do not 
employ the weak boson fusion cuts, opting instead for a selection criteria that give sizable 
cross-sections for pp — > W + W~jj. Our choice of cuts is inspired by those that are made 



in the first analyses of tt production by ATLAS and CMS collaborations [6C, 61]. We do, 
however, plot distributions which are interesting in the context of reducing the W + W~jj 
background to the Higgs-boson searches in weak boson fusion. For example, we study the 
relative jet rapidity At^v, = r\j x — rjj 2 distribution and the opening azimuthal angle of the 
two leptons ^ 1 ; 2 . Given that the center-of-mass energy of collisions at the LHC after the 
longer shutdown at the end of 2012 is not fully decided yet, we also find it interesting to 
show the behaviour of the cross-section as a function of y/s. 

We consider proton-proton scattering pp — > W + W~ jj at center-of-mass energy ^fs = 
7 TeV. We impose the following cuts, inspired by tt searches at the LHC: 
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Figure 9: In the left pane, we show the production cross-section of the process pp — > (W + — > 
v^ii + ) [W~ — > e~v e )jj at the 7 TeV run of the LHC in dependence of the factorization and 
renormalization scales fiF — M-R — Mi a t LO and NLO in perturbative QCD. In the right pane, the 
dependence of the cross-section on the center-of-mass energy i/s is shown. LO results are shown in 
dashed blue; NLO results are in solid red. Three choices of \i are shown: [i = Mw, 2Mw, 41%- 



jets are denned using the anti-A:^ algorithm [g^] as implemented in Fast Jet [63], with 
AR jlj2 = yj (77ji - r] 32 f + (^-i - 4> l2 ) 2 > 0.4; (3.1) 



• jets are required to have transverse momentum p±j > 30 GeV and the rapidity 
M <3.2; 

• charged leptons are required to have transverse momenta > 20 GeV and the 
rapidity \r]i\ < 2.4; 

• missing transverse momentum is required to satisfy p_i_ jm i ss > 30 GeV. 

In the left pane of Fig. ^ we show the dependence of the cross-section pp — > W + W~ — > 
[i^v^eVf, jj at the 7 TeV run of the LHC, on the factorization and renormalization scales, 
which we set equal to each other. At leading-order, the cross-section falls with the scale 
fi, which is attributable to the behaviour of the strong coupling a s . Considering a range 
of factorization/renormalization scales M\y < /i < 4M^y and choosing the central value 
fi = 2M\y, we obtain a cross-section cjlo = 46 ± 13 fb. At next-to- leading order, the 
dependence on /x is dramatically reduced and the cross-section becomes onlo = 42 ± 1 fb. 
Such a decrease in the scale dependence is typical of NLO results, and indeed one of the 
primary motivations for performing calculations at next-to-leading order in pQCD. At the 
scale [i = 2Mw, the NLO corrections increase the cross-section by about 2%. Assuming 
fifty percent efficiency, with 5 fb -1 of data at the 7 TeV run of the LHC, we expect about 
400 dilepton events e + /i~, e/i + , e + e~, fi + ^ ■ 

It is interesting to know how the cross-section for W + W~jj production changes with 
the collision energy. In the right pane of Fig. 0, we show that the dependence of the 
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Figure 10: Kinematic distributions for jets in the process pp — > u fl p + e~P e jj at the 7 TeV run of 
the LHC at LO and NLO in perturbative QCD. The bands show uncertainty on the renormalization 
and factorization scale p, for Mw < p < 4Mvy, while the lines show results for p — 2Mw- 



NLO cross-section on the center-of-mass energy ^fs is very close to linear. Again, the 
significant reduction in uncertainty in the NLO prediction for the cross-section is obvious 
from Fig. ^. It follows from Fig. [|| that the optimal 2 renormalization/factorization scale 
increases with the center-of-mass energy smoothly interpolating between p = 2Myy at 
7 TeV and p = 4 M w at 14 TeV. 



We now turn to the discussion of kinematic distributions. In Fig. [10| we show the trans- 
verse momentum distribution of the hardest and next-to-hardest jets and the distribution 
of the total transverse energy Hj^tot- For all distributions, the scale dependencies are 
reduced and shapes of the distributions are, typically, not distorted. Note, however, that 
the NLO QCD corrections make the jet transverse momenta distributions and the i?r,TOT 
distributions somewhat softer, which is caused, at least partially, by our use of a constant, 



2 We define the "optimal" renormalization/factorization scale as the value of p for which next-to-leading 
order corrections are the smallest. 
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Figure 11: Kinematic distributions for leptons in the process pp — > v^^ + e"v e jj at the 7 TeV run of 
the LHC at LO and NLO in perturbative QCD. The bands show uncertainty on the renormalization 
and factorization scale /x, for My/ < A* < 4Mvy, while the lines show results for [i — 2Mw- 



rather than a dynamic, renormalization scale in the LO calculation. We show lepton kine- 
matic distributions in Fig. 11. Similar to jet distributions, lepton transverse momentum 
and the missing energy distributions are softened by the NLO QCD corrections. 

A few other distributions which are relevant for designing cuts for Higgs searches are 
presented in Fig. The distribution of the relative azimuthal angle between the two 
leptons is peaked at (f> e -^+ = vr, with the NLO corrections making almost no change in 
the shape of that distribution. The pseudorapidity difference between two leading jets, 
defined as Arjjij2 = Vji — ilj2, peaks at small values of Arjjij2 and falls off rapidly for 
larger values. The invariant mass of the leptons and the transverse mass of the l^-bosons 3 
become somewhat softer once the NLO QCD corrections are included. A discussion of how 



these distributions can be used in searches for the Higgs boson can be found in Refs. pq, 59 



64]. The availability of NLO QCD predictions for those distributions should, potentially, 
improve the reliability of such analyses since, as follows from the discussion in this paper, 
theoretical uncertainties are reduced considerably. 



4. Conclusions 

In this paper, we computed the NLO QCD corrections to the production of a W + W~ pair in 
association with two jets in hadron collisions. We only considered the QCD contribution to 
this process, ignoring the possibility that it can also occur through exchanges of electroweak 
gauge bosons. Our calculation includes the leptonic decays of VF-bosons and accounts for 
all spin correlations exactly. 

3 We define Mj_ iWW = (E Xtl+l - +£±, miss ) 2 - (p_i,i+i- -P±,mi ss ) 2 , where £j_, miss = Jp± imiBB + mf +l _. 



-17- 





50 100 150 200 250 300 
m e ^ + [GeV] 



0.00 



100 200 300 400 500 600 700 
M ± , ww [GeV] 



Figure 12: Distributions of jet pseudorapidity difference, lepton opening angle and invariant 
masses for pp —> v^^i + e~v e jj at the 7 TeV run of the LHC. LO results are shown in blue; NLO 
results are in red and orange. The uncertainty bands are for scale Mw < H < 4Mw, and the solid 
lines show the results at /i = 2Mw- 



The computation of NLO QCD corrections was performed using the method of D- 
dimensional generalized unitarity & [|. Practical implementations of the generalized uni- 
tarity technique require color ordering 4 ; for this reason, the presence of any colorless particle 
leads to additional complication since colorless particles can not be ordered. Most processes 
for which the NLO QCD corrections have been computed using the on-shell methods in- 
volve at most one colorless particle. The results of this paper and of Ref. |24[] show that 
generalized unitarity methods can be efficiently used to deal with processes with a larger 
number of colorless particles, although the most general framework for that is yet to be 
understood. 

We studied some phenomenology of the W + W~jj production at the Tevatron and 
4 See, however, a recent discussion in Ref. . 
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the LHC, using y/s = 7 TeV for the center-of-mass collision energy of the latter. We also 
explored the behaviour of the NLO QCD cross-section for pp — > W + W~jj as a function 
of the center-of-mass energy at the LHC and find that, to a good approximation, the NLO 
cross-section grows linearly with the energy of the collider. For the renormalization and 
factorization scales set to \i = M\v and \i = 2M\y at the Tevatron and the LHC respectively, 
the radiative corrections for both colliders are moderate; in fact they are very small for 
collisions at 7 TeV. We show that the uncertainty in the theoretical prediction, estimated 
by changing factorization and renormalization scales in the range 0.5Myy (Myy) < fi < 
2M W (4M W ) at the Tevatron (LHC) is better than 10% if the NLO QCD corrections 
are included. Of course, at that level of precision other uncertainties - such as e.g. the 
imperfect knowledge of parton distribution functions - become important. We considered a 
number of kinematic distributions that involve lepton and jet momenta and observed that 
energy-related distributions (p±, H^ot) become softer once the NLO QCD corrections are 
included and that shapes of angular distributions are hardly affected. We also discussed the 
significance of pp — > W + W~jj process as an irreducible background for the production of 
the Higgs boson in association with two jets at the Tevatron, as well as kinematic variables 
useful to disentangle a Higgs signal from the W + W~ background. 
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A. Results at a fixed phase space point 



In this Appendix, we shall give numerical results for some of the tree-level, primitive and 
full virtual amplitudes used in this calculation. For the sake of brevity, amplitudes for 
some flavors and helicities are not been reported here. However, we also give results for 
squared amplitudes, summed over helicities and color. 

We begin by considering the process — > (qq) + (W + — > + fi + ) + (W~ — > e~ + D e ) + 
g + g, and use the phase space point defined by the following values of momenta 

pf = (-500.00000000000000, -500.00000000000000, 0.00000000000000, 0.00000000000000), 
Vl = (-500.00000000000000, 500.00000000000000, 0.00000000000000, 0.00000000000000), 
pi" = (85.5312248384887, -8.22193223977868, 36.16378376820329, -77.0725048002413), 

= (181.42881161004266, -57.85998294819373, -171.863734086635, -5.611858984813), 
pf = (82.84930107743558, -65.90954762358915, -49.89521571962871, 5.51413360058664), 
pf = (381.47038530081545, 190.18527704151887, 292.042940984587, -155.113300136598), 
p 9 7 = (54.23140701179994, -31.13301620817981, -7.9279665679114, 43.69128236111634), 
p% = (214.48887016141776, -27.06079802177751, -98.519808378615, 188.59224795994947). 

Our convention for displaying four-momenta is p = (E,p x ,p y ,p z ); all momenta are given 
in GeV. 

We only include results for the case in which the helicities of the uu are (+, — ), even 
though the opposite helicities do contribute via an intermediate vector boson. We also 
do not include the results for dd - these can be obtained by switching the order of the 
VF-bosons and modifying the 7/Z couplings in equation (2.4). In Tables ||, ||, || and we 



(A.l) 



give tree-level amplitudes as well as the ratios of the unrenormalized virtual amplitudes to 
the tree-level amplitudes 

1 Ai [1/2] 1 4 V2] ( A 2) 

r x = — r [ / = ij — , (A.2j 

c r A c r A 

where cr = S -=7- 1, and the renormalization scale is iir = 150 GeV. The tree- 

(47r) /-< T(l — 2e) 



level amplitudes Aq are defined in Eq. (2.1), while the primitive amplitudes A\ are defined 



in Eq. (2.7). The one- loop amplitudes are calculated in the four-dimensional helicity scheme 



1 66 , |67| . Finally, in Table [| we give the ratio 

„ 4yr E {hc i } Re(A t -A 1L *) 

bA = 1 yi tree 1 2 ' ( A - S > 

a s E{hel} \A tTce \ 2 

where the sum is over all helicities for the quarks and gluons. 

We now consider the case of — > (qq) + {W + — > v ^ + p + ) + (W~ — > e~ + i7 e ) + q 3 + q±. 
We use the same momenta as in equation ( |A.l| ), with the modification that the last two 



momenta in Eq. ( A.l ) are now those of a q^q^ pair, p^ — > p^ 3 , p| — > pg 4 . For the sake of 
brevity, we restrict the results given here to two sets of flavors: uucc and udsc (we are thus 
working with an "s-amplitude" ) . The flavor structure of the first set is given in Eq. ( 2,1C| ) . 
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Amplitude 


1/6 2 


1/6 


e° 


^(9i + )03~jfl4>?2 ) 

r i{q~i + ,g% ,gl ,q2~) 


-1.000000 


2.294240 - i 3.141593 


Q Q/M 1 l n ni OOH7 

— o.o441od + z y.yizzu/ 
0.4601166 + ^2.774496 


MVl\ ,9s ,94 ><?2 ) 
ri(qi + ,9s,gt,q2) 


-1.000000 


2.294240 -i 3.141593 


U. ( UODoll -p 2 D.DOZD4U 

0.3739239 + ^2.687541 


A {qi + ,gti94,<i2) 

ri(q~i + ,g£,g4,q2) 


-1.000000 


2.294240 -i 3.141593 


-5.998084 -i 5.572010 
0.5484790 + i 3.010535 


A (qi + ,gt,gt^2) 
n(qi + ,g£,gi,q2) 


-1.000000 


2.294240 - i 3.141593 


-10.07279 - i 3.926576 
0.4741562 + ^2.846111 



Table 1: Numerical results for the primitive tree- level amplitude ^4o(<Zi, 53, <?4, 92), in units of 
1CP 10 GeV 4 and the ratios of primitive one-loop amplitudes r\(q~\, 53, 54, 92)- 



Amplitude 


l/e 2 


1/6 


6° 


A o{q~i + ,g% ,q2 ,gl ) 
r i(qi + ,g3,q2,g4~) 


-2.000000 


2.993440 + i 0.000000 


-5.097350 + i 3.386328 
-0.07739397 + i 3.420824 


A (q~i + ,g3,q2,gt) 
n{.qi + ,g^ ,q% ,g£ ) 


-2.000000 


2.993440 + i 0.000000 


-4.426865 + i 4.803504 
6.347479 + i 5.196425 


A (qi + \g£ \%. ',9a) 
ri{q~i + ,gt ,q2 ,g~4~) 


-2.000000 


2.993440 + i 0.000000 


-4.749089 + i 1.306764 
-0.8538774 + ^3.373345 


Ao(<Ti + ,gt ■>% i9t) 

ri(q~i + ,gt ,q2 ,gt) 


-2.000000 


2.993440 + i 0.000000 


-8.206743 + i 2.583236 
6.051784 + i 4.612948 



Table 2: Numerical results for the primitive tree-level amplitude Ao(gi, 53, q%,g4), in units of 
10 -10 GeV -4 and the ratios of primitive one-loop amplitudes ri(q\, 53, q2, 54)- 



We shall also restrict ourselves to the helicities qqqq = (+,—,+,—), although for the former 
set of flavors, there are four different helicity combinations that are used in the calculation. 
We give the ratios 



1 B 



W 
1 



p (A.4) 

for i = a,b,c,d, where Bq is defined in Eq. and b[^ are defined in equations ( 2,12j ), 
( |2.13|) , and The results are shown in Tables |6| and 0. We also give the ratios 

q 4vr E {he i } Re(^-g 1L -) 

bB = IKtrccl2 ^ A - 5 ^ 

a s E{hol} \ B I 

in Table § 
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Amplitude 


l/e 2 




e° 


A ( V — — — \ 

r i{qi + ,q2 ,93 ,9a) 


-3.000000 


-0.9503441 -i 3.141593 


8.44153b — i 13.29854 
-6.047837 - i 9.654414 


Ao(qi + ,q 2 ,93,gi,) 

n(qi + ,q 2 ,93,gi) 


-3.000000 


-0.9503441 -i 3.141593 


o. ( 21oo4 — I 11.48614 
0.9335325 - i 8.464906 


Ao(qi + ,q2,9t,9A) 
n(qi + ,q2,93,9A~) 


-3.000000 


-0.9503441 -i 3.141593 


10.74717 + i 4.265245, 
-6.036407 -i 10.58605 


Ao{qi + ,q2,93,9t) 

ri{qi + ,q 2 9t,9A) 


-3.000000 


-0.9503441 -i 3.141593 


18.27953 + i 1.343340 
0.3979266 -^9.181091 



Table 3: Numerical results for the primitive tree-level amplitude Ao(q\, 92,53,54), in units of 
10~ 10 GeV~ 4 and the ratios of primitive one-loop amplitudes ri(q\, q 2 , <?3, gi)- 



Amplitude 


e° 


M(qi + ,q 2 ,93,gA ) 

[1/2], _ + - - — \ 

n + ,q 2 ,9 3 ,9a ) 


8.441536 - i 13.29854 
(-0.3523178 - i 4.071390) x 10~ 2 


Ao(qi + ,q 2 ,93,9t,) 
r [ i /2 \qi + ,q2,93,gi) 


3.721334 - i 11.48614 
0.000000 + i 0.000000 


A (qi + ,q 2 ,gi,g4) 

[1/2]/ _ + _ + _x 

n {qi + ,q 2 ,93,9a ) 


10.74717 + i 4.265245 
0.000000 + i 0.000000 


Ao{qi + ,q 2 ,gt,9t) 

[1/2]/-+ - + +x 

n (qi + ,q 2 93,9A) 


18.27953 + i 1.343340 
(-3.142652 + t 1.567695) x 10~ 2 



Table 4: Numerical results for the primitive treedevel amplitude Ao(qi, q 2 , 53, 54), in units of 
1CP 10 GcV~ 4 and the ratios of primitive onedoop amplitudes r^ 2 ' (gi, q 2 , 53, 54). There arc no 
singular contributions from these onedoop amplitudes, and swapping the gluons simply changes the 
sign of the amplitude. 
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Ratio 


iA 2 


iA 


e° 


Z\A tree (d,d,g,g)\ 2 






9.887737 x 10~ 20 


S A {d,d,g,g) 


-8.666667 


-2.836720 


-0.6913131 


£|-4 trcc (^,5,2)| 2 






3.743231 x 10" 20 


SA(u,u,g,g) 


-8.666667 


-2.786885 


-4.673601 



Table 5: Numerical results for the tree-level amplitude squared, in units of GeV 8 , and the ratio 
of virtual over tree-level squared amplitudes Sa summed over all helicities and colors. 



Amplitude 


iA 2 


iA 


e° 


B (u,u,c, c) 
r a (u,u,c, c) 
r a (u,u,c,c) 
r b (u,u,c, c) 
r c (u,u,c, c) 
r d (u,u,c, c) 


-2.000000 
-2.000000 
-1.000000 
-1.000000 


3.066474 + i 0.000000 
4.119961 + i 0.000000 
2.294240 -i 3. 141593 
-3.350152 - i 3.141593 
-0.6666667 + i 0.000000 


0.6391654 + i 5.544406 
2.658086 + i 2.684586 
3.634715 + i 2.090514 
0.1918562 + ^2.854994 
-3.028899 - i 10.77523 
-2.301323 - i 1.838568 



Table 6: Numerical results for the primitive tree- level amplitude B (u,u,c,c), in units of 
10~ 10 GcV~ 4 and the ratios of primitive one-loop amplitudes r». 



Amplitude 


l/e 2 


iA 


e° 


B (u,d, s, c) 
r a (u,d, s, c) 
r a (u,d, c, s) 
r b (u,d, s, c) 
r c (u, d, s, c) 
r d (u,d, s, c) 


-2.000000 
-2.000000 
-1.000000 
-1.000000 


3.066474 + i 0.000000 
4.119961 +i 0.000000 
2.294240 -i 3.141593 
-3.350152 -i 3.141593 
-0.6666667 + i 0.000000 


0.3350897 -i 0.6484033 
-7.426922 - i 0.3913681 
-7.135027 - i 13.92234 
0.7221019 +i 6.182924 
-7.635220 - U0.639296 
0.9665113 - i 2.094395 



Table 7: Numerical results for the primitive tree-level amplitude B (u,d, s,c), in units of 
10~ 10 GcV~ 4 and the ratios of primitive one-loop amplitudes rj. 
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Ratio 


1/6 2 


1/6 


e° 


£|£ trce (u,u,c, c)| 2 






1.037139 x 10~ 21 


S B (u,u,c, c) 


-5.333333 


7.587051 


5.395242 


Z\B tTee (u,d,s,c)\ 2 






1.123763 x 10~ 23 


S B (u,d,s,c) 


-5.333333 


7.587051 


-15.91575 



Table 8: Numerical results for tree-level amplitudes squared, in units of GeV 8 , and the ratio of 
virtual over tree-level squared amplitudes Sb summed over all helicities and colors. 
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